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A theory for the photon statistics of a random laser is presented. Noise is described by Langevin 
terms, where fluctuations of both the electromagnetic field and of the medium are included. The 
theory is valid for all lasers with small outcoupling when the laser cavity is large compared to the 
wavelength of the radiation. The theory is applied to a chaotic laser cavity with a small opening. 
It is known that a large number of modes can be above threshold simultaneously in such a cavity. 
It is shown the amount of fluctuations is increased above the Poissonian value by an amount that 
depends on the number of modes above threshold. 
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I. INTRODUCTION 

A random laser is a laser where the necessary feedback 
is not due to mirrors at the ends of the laser but due to 
random scattering inside the medium jl], |], It was 
long argued how to distinguish such a random laser from 
a random medium with amplified spontaneous emission 
(ASE) — in the former, the randomness is essential for 
providing feedback, whereas in the latter, scattering only 
increases the dwell-time in the medium and thus the am- 
plification factor. Two years ago, the first experimental 
proof of a random laser was given Q. It was demon- 
strated that the lasing action was indeed due to the ran- 
domness of the medium, by measuring the emitted radi- 
ation at different points on the surface of the sample and 
showing that the peaks in the radiation spectrum were 
completely different at different points. 

Earlier experiments |5|, |^, Q were only able to prove 
ASE in random media, frequently referred to as "laser- 
like emission" . In a medium with saturation both laser 
action and ASE lead to a dramatic narrowing of the emit- 
ted light profile upon crossing some threshold so that this 
criterion does not necessarily signal a laser. Most "tra- 
ditional lasers" are characterized by emitting coherent 
radiation above threshold so that considering only the 
intensities and forgetting about the fluctuation proper- 
ties is insufficient. Recently the first two measurements 
on the photon statistics of a random laser have been pub- 
lished. The group of Papazoglou reports that the emit- 
ted radiation becomes only partially coherent j^] whereas 
the group of Cao reports that the statistics become com- 
pletely Poissonian M. 

The theoretical description of random lasers has in the 
past focused on the light intensity inside in the laser. 
Photons were considered as classical particles that diffuse 
or move in some other way repeatedly through the sample 
while being amplified. (The literature on this and simi- 
lar methods is numerous, some more general, some more 
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focusing towards a particular system; see e. g. Ref. 
for one of the earlier papers.) In this way the intensity 
of the emitted radiation can be computed, confirming 
the observed narrowing of the emission line far above 
threshold. No results for the fluctuation properties, how- 
ever, can be derived in this way. Recently, random lasers 
are also simulated by the finite-difference time domain 
method (FDTD) pJj . While this method in principle 
can incorporate quantum fluctuations on a microscopic 
level, the computational effort is prohibitively large, so 
that at most two-dimensional samples can be treated (see 
e. g. Ref. |l2|), and most of its value is for one-dimensional 
applications (see e.g. Ref. 13). Furthermore, only short 
time series can be computed with acceptable effort so 
that the fluctuation properties of the emitted radiation 
are not accessible. A different, analytical, approach to 
noise in random lasers has recently been put forward by 
Hackenbroich et al. |Q . Since they do not include mode 
competition, their work is only applicable near threshold. 

For a linear medium, i. e. a medium where, in contrast 
to a laser, saturation effects can be neglected, the statis- 
tics of the emitted radiation can be computed directly, 
e.g. by the method of input-output relations jig] . No 
theory of comparable power exists for lasers. The theo- 
retical treatment of "non-trivial" lasers has in the past 
focused on the Pctcrmann factor (see Refs. |l6|, 17, [l8| 
for a definition). It is a geometry- related factor that 
describes by how much the excess noise of the emitted 
radiation is larger than for a "simple" single-mode laser 
- assuming that the non-trivial laser behaves the same 
way as a single-mode laser, which is basically equiva- 
lent to neglecting mode competition effects. (It should 
be stressed that the Petermann factor only gives infor- 
mation about the radiation far above threshold; it gives 
no information on threshold behavior.) Since the Peter- 
mann factor is a geometrical factor it can be computed 
for a linear medium and then used for the corresponding 
system filled with a medium with saturation. The Pe- 
termann factor has been derived for arbitrary geometries 
(see e.g. Ref. [l9[ po|) but also random media could be 
treated Hlgf]. 

There thus is a need for a theory that allows one to 
compute the photon statistics of the emitted light for 
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"non-trivial" lasers, in particular this includes random 
lasers. In this paper such a theory based on Langevin 
terms, also referred to as Langevin noise sources, is pre- 
sented. Langevin terms have successfully been used to 
describe the radiation properties of linear media from a 
microscopic model |2fj| . On a higher level, they were 
used to describe random linear amplifying media j24j] 
where the Langevin terms included both fluctuations of 
the electromagnetic field and sample-to-sample fluctua- 
tions of the properties of the random medium. None of 
these theories included saturation effects of the medium 
so that they break down when the lasing threshold is 
approached. Apart from saturation effects for a single 
mode, a large number of modes can be above threshold 
simultaneously so that mode-competition is impor- 
tant and cannot be neglected. 

This paper is organized as follows: In Sec. II the model 
for the photon statistics inside the laser is described and 
the model equations are derived. These are then solved in 
Sees, g and Sec. adds the necessary modifications 
to go from the fluctuations inside the laser to the fluc- 
tuations of the photocurrent emitted by the laser. Until 
this point all results are valid for arbitrary lasers, pro- 
vided that the outcoupling is weak and the volume of the 
lasing medium is much larger than the cube of the wave- 
length. In Sec. VI we show how to apply the formalism 



developed in this paper to three exemplary systems and 
demonstrate thereby that it can indeed desc ribe all rele- 
vant properties of lasing action. In Sec. VII the random 
laser is treat ed, and its photon statistics are computed. 
In Sec. VIII we try to explain the experimental results 
mentioned above. We conclude in Sec. [X. 
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II. MODEL 



FIG. 2: Amplification is modeled by a four-level system, 
where lasing action (marked by the wiggled line) is from the 
third to the second level. Dashed lines mark transitions that 
are much faster than the other ones and thus need not be 
included in the description. 



medium is modeled by a four-level laser dye (see Fig. ||), 
where the lasing transition is from the third to the second 
level. The transition from the second level to the ground 
level is assumed to be so fast that the second level is al- 
ways empty. The density of excited atoms (i. e. atoms in 
the third level) at point r in the cavity is N(r). Excita- 
tions are created by pumping with rate P(f) and can be 
lost non-radiatively with rate a(f). 

Coupling between the electromagnetic field and the 
medium depends on two quantities, namely the eigen- 
mode profile 0j(r) of mode i, and the transition ma- 
trix element w{u>) of the atomic transition 3^2. [Fre- 
quently w(lo) will be a Lorentzian centered around some 
frequency Q.] The coupling of mode i to the medium at 
point r is then given by Ki{r) = w(uji)\Qi(r)\ 2 . 

The semiclassical equations of motion for rti and N(r) 
are (the time argument for all quantities has been sup- 
pressed) 



We consider a optical cavity that is coupled to the 
outside by an opening that is small compared to the 
wavelength of the emitted radiation (see Fig. [l]). Since 
the opening is small, there exist well-defined modes in 
the cavity, each with a well-defined eigenfrcquency Wj, 
i = 1,... , N p , and a eigenmode profile ®i(r), and all 
modes are non-overlapping |3^| . (In the language of ran- 
dom lasers, this is a "resonant-feedback laser".) Each 
mode i thus can be described by the number m of pho- 
tons in it. Photons in mode i can escape through the 
opening with rate gi. 

The cavity is filled with an amplifying medium. The 




FIG. 1: A (chaotic) cavity is coupled to the outside via a small 
opening. The cavity is filled with an amplifying medium. The 
light emitted through the opening is detected. 



hi = -g l n i + J d 3 r (m + l)Ki(r)N(r) , (la) 

N(f) = P(f) - a (r)N(r) - ^](n t + l)Ki(f)N(f) . 

i=i 

(lb) 

"Semiclassical" means that all emission events, pumping 
events, . . . are assumed to be deterministic, with spon- 
taneous emission described by the addition of a virtual 
photon to rii when computing the transition rates fl35| . 

To include the randomness of all processes, Langevin 
terms have to be added to Eq. |f]). The four random 
processes are the escape of photons (described by the 
Langevin term r,-), pumping [described by 3>(r)], relax- 
ation of the medium [described by a(r)}, and emission of 
a photon into mode i at point f [described by ^(r)]. 
Each of these terms has zero mean, and a correlator 
that follows from the assumption that the elementary 
stochastic processes have independent Poisson distribu- 



3 



tions, hence 



{r i (t)T j (t')) = 8 ij S(t-t')g i (n i ) , 



(2a) 



{a{r,t)a(r',t')) = 8 3 (r - r')8{t - t')a(r)(N(r)) , 

(2b) 

($(r, t)*(r t')) = 8 3 (r- r')8(t - t')(P{r)) , (2c) 

(^(r.t^f ',*')) = %^ 3 (r - r')<5(t - t') x 

i^(r)((n t + l)7V(r)) . (2d) 



Eq. ( 2d ) corresponds with the correlator given in Eq. (5b) 
of ReT||. 

Adding the terms from Eq. (^J) to Eq. (|l|) gives the 
complete equations of motion: 



hi = -g ini +T t + J d 3 r(n z + l)Ki(r)N(r) 

+ J d 3 rVi(r) (3a) 
N(r) = P(r) + $(f) - a(r)N(f) + a(f) 

N p N p 

-J^(n t + l)K t (r)N(r)~J^^ l (f) (3b) 



The sign of the Langevin terms may be chosen freely 
as long as the term \E'i(r) has the opposite sign in the 
equations for h and N. 



average solution: 

Shi = -gMi + J d 3 f{hi + l)Ki(r)6N(r) 

d z r5n i K l {rjN(r) + J d 3 r^i{r) (4a) 
SN(f) = $(r) - a(r)8N(r) + a(f) - ^(n 4 + l)Ki(r)SN(r) 

i 

_ _ Np 
- SniKi (f)N(f) -J2*i (*o ( 4b ) 

i i—1 

For convenience, we will label the sum of the Langevin 
terms in Eq. (pia| ) as fi, and the sum in Eq. (|l]) as g(r). 
Evaluating the Langevin terms from Eq. (^|) at the aver- 
age solutions h and TV gives 



(fifj 



gi ni+ / d 3 r{h t + l)Ki{r)N{r) 



= ^SijgiHi 



(5a) 



(g(r)9(r')) = 8 3 {r- r ') a{r)P(r) + N(r) 



Np 



.^(fti + l)K t (r)N(r) 



28 3 {r-r')P(r) 



(M^) = -(h l + l)K l (f)N(f) 



(5b) 
(5c) 



III. LINEARIZATION 

Eq. (|^) cannot be solved by direct numerical meth- 
ods since Langevin terms cannot be represented as "real" 
numbers. The only practicable way to proceed is to lin- 
earize the equations. First, we write rii = hi + Srii and 
N(r) = N(r) + SN(r), where hi = (n») and N(f f ) = 
(N(r)) are the average solutions. We assume that these 
average solutions are identical to the solutions of the de- 
terministic rate equation (0). This is equivalent to the 
factorizing approximation (mN(f)) f=a (n,) (N (r)) . For a 
single-mode cavity like used in cavity-QED this is a bad 
approximation, leading to errors of up to a factor 1/4 in 
the computed average photon density, but if the number 
of modes in the cavity is large — which is the case that 
we are interested in — this factorization is valid |28| . 

Inserting this solution, Eq. ([|) can be reformulated so 
that only drii and 8N(r) remain as variables. Lineariza- 
tion means that only terms proportional to dtii or 8N(r) 
are kept, i. e. terms proportional to 8rii8N(f) are omit- 
ted. (This is justified as long as the variance is sufficiently 
smaller than the mean. This condition is equivalent to 
the validity of the factorizing approximation used above. 
It can be checked self-consistently from the computed 
results.) This way one arrives at an equation for the 
fluctuations alone, where the coefficients depend on the 



IV. DISCRETIZATION AND NUMERICAL 
SOLUTION 

We now discretize the equations in space by picking 
points fj, j — 1,... ,N S . Defining Kij = Ki(fj) and 
Nj = N(fj) (analogously for all other quantities), the 
stationary densities hi and Nj from Eq. (Q) are the so- 
lution of the equations 

gthi = Y^(hi + ^KijNj (i = 1, . . . , N v ) (6a) 



i=i 



N s 



Pi = ajNj + ^(hi + l)K l3 N, (j = 1, . . . , N s ) 

1=1 

(6b) 

This equation cannot be solved analytically but a nu- 
merical solution is straightforward (even though it may 
be numerically expensive if N p and/or N s are large). 
Eq. (^) now becomes a linear ODE, 

d / Sn t 

-ft + T.jKijNj (hi + l)K i} 

-KijNj -a, - EiC^i + l)K iS J V 6N i 



9j 



(7) 
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where it is understood that all indices i run from 1 to 
N p and all indices j from 1 to N s , so that the previous 
equation can be written as a (jV p + iV s )-dimensional ma- 
trix equation 5fif = ASM + C. Computing from A its 
matrix U of eigenvectors and its vector £ of eigenvalues, 
the formal solution can immediately be written down as 

i...n p +n b t 

8Nj{t)= Y, / dt'Ujke^-^U^Ciit'). (8) 

k,l 

Since the vector C consists of Langevin terms, a numer- 
ically computed solution of Eq. m) is not meaningful. 
Instead of 6Afj(t) alone one has to consider correlators 
(5Afj(t)SAfji (t)) . Noting that the C's are delta-correlated 
in time, and we are interested in t — > oo (as we are not in- 
terested in intermittent behavior when switching on the 
laser) we arrive at 



-cT/2 



1-N P +N JJ-llJ-l 

mwr) = - E J t — wd 



h ! in n 



(9) 



Inserting the expectation values of the correlators from 
Eq. (|^) gives the final result where a numerical solution is 
easy once the average solution n,, Nj is known. ((£&£;) 
has to evaluated at the average solution and does thus 
not depend on time.) 



V. OUTCOUPLING 

So far we have considered the number of photons rii in 
the i-th mode inside the cavity. For practical purposes 
one is more interested in the photo current I emitted from 
the cavity. (J gives the number of photons emitted per 
unit time and is thus equal to the photon flux integrated 
over the entire cross-sectional area.) Even though the 
photons from different modes i are emitted through the 
same opening, each mode has a distinct frequency uii so 
that the modes are easily distinguished on the outside. 
We can thus define the photo current ji (t) = ji + Sji (t) 
through the opening due to the i-th mode in the cavity. 
The photo current can for example be measured by an 
(ideal) photodetector that absorbs the emitted photons. 
The fluctuations of the photo current within some time r 
(we assume the limit r — * oo) are quantified by the noise 
power 



P 



1 
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lim - / dtdj(0)Sj(t) 
t^oo r J_ T / 2 



(10) 



The ratio Ti = Pi/ji is called the Fano factor and is 
frequently used to describe the fluctuation properties of 
optical radiation. 

In Sec. II we have introduced the loss rates gi. From 
their definition it is obvious that the mean photo current 
ji is 



9irii 



(11) 



FIG. 3: The loss rate of photons inside the cavity is given by 
the ratio of the probability ti that a photon incident on the 
outcoupling mirror is transmitted and the time T needed for 
one round trip through the cavity. The photons emitted from 
the cavity are detected by an ideal photodetector. 



To also compute the fluctuations Sji we need to treat 
the outcoupling in more detail. In a traditional laser 
(see Fig. ||) the loss rate is given by the ratio of the 
transmission probability ti (in classical optics referred to 
as "transmittivity" ) through the outcoupling mirror and 
the round-trip time T through the cavity, 



gi = U/T 



(12) 



The transmission through the outcoupling mirror 
changes the noise of the signal compared to the noise 
inside the cavity, and the Fano factor of the emitted ra- 
diation is fE6| 



(13) 



This equation can, apart from following the quantum- 
optical approach of Ref. [2(| also be understood by the 
following simple argument: The fraction (5rii5rii)/fii on 
the right-hand side is the Fano factor of the radiation 
trapped inside the cavity in mode i. With probability 
ti the detector will "see" the radiation inside the cavity, 
and with probability 1 — ti it will sec reflected vacuum 
fluctuations (which have a Fano factor equal to 1). 

The Fano factor for a measurement where the photons 
emitted from the cavity in all modes are detected simul- 
taneously is 



T = 



(14) 



It is immediately obvious that ti and gi can for a tra- 
ditional laser be identified by properly choosing the unit 
of time (for the simple laser from Fig. |j: b y ch oosing T 
as the unit of time). We will show in Sec. |VI]| that this 
is also possible for a random laser. In the following when 
giving numerical values or distribution functions for 
this identification has been made. 



VI. COMPARISON OF LASING REGIMES 

To demonstrate the application of the formalism pre- 
sented in this paper and the validity of the approxima- 
tions made in this paper we first want to discuss three 
simple cases not involving random media. For simplicity 
we set a = w = 1, N s = N p and K = const. This reduces 
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FIG. 4: Comparison of the Fano factor T for three different conditions. The left axis (solid line) depicts the Fano factor for the 
integrated emitted radiation, the right axis (dashed line) for the lasing mode only, (a) Laser with just a single mode, (b) Laser 
with a cavity supporting 10 modes where one mode is coupled out much less than the others, thus effectively modeling a 
single-mode laser with /3 « 0.1. (c) Laser with 10 identical modes. 



the number of parameters significantly without reducing 
the physical content. 

The physical features of a laser (in contrast to a linear 
amplifier) are easily understood in the following picture: 
A certain number of excited atoms are created by pump- 
ing within a certain time, and each of those excitations 
has to be "consumed" either by nonradiative relaxation 
or by emitting one photon from the cavity. For high pho- 
ton number in the cavity, nonradiative relaxation can be 
neglected, and each pumping event eventually leads to 
the emission of one photon from the cavity. The fluc- 
tuations of the integrated photo current are thus equal 
to the fluctuations of the pump source, assumed to be 
Poissonian throughout this paper. 

In Fig. ||(a) the single-mode laser (N p = 1) with a small 
opening (g = 10~ 2 ) is treated. The computed curve re- 
produces the features of a "traditional" laser. The precise 
location of the maximum is somewhat off (see the discus- 
sion of the factorization approximation above, or refer 
to Ref. [^7j for a more detailed discussion of the effects 
of different approximations on the computed curve near 
the lasing threshold) but its height reproduces the exact 
quantum-mechanical value well. For high values of the 
pumping, the photon statistics of the emitted radiation 
becomes Poissonian, as qualitatively explained above. 

In Fig. |J(b) we have modeled a laser with one mode 
coupled to the outside with g = 10~ 2 and the other 
9 modes with g = 1CP 1 , hence iV p = 10. (The value 
g = 10~ 2 was chosen for scaling the axes of the figure.) 
The mode with the smallest g will be the lasing mode, 
whereas radiation in the other modes quickly escapes 
to the outside so that no significant number of photons 
can accumulate in those modes. This basically models a 
single-mode laser where only a fraction (3=1 /N p of the 
spontaneous radiation is emitted into the lasing mode. 
((3 is called the spontaneous emission factor. An ideal 
cavity-QED laser has (3 = 1 whereas a semiconductor 
laser can have a 8 as low as (3 = 10~ 8 .) The behavior 
is similar to Fig. [|(a), except that the peak of the Fano 
factor of the lasing mode is larger by about a factor 8. 
For small beta, one expects a scaling cx (3~ 1 / 2 w 3 pgj ] 
but (3=1 and (3 = 0.1 are too large for that scaling to 



be exactly valid. 

In Fig. ^(c) the system is kept at N p = 10 with all 
gi = 10~ 2 . The total radiation depicts the same qualita- 
tive behavior as for the two cases presented so far but the 
radiation emitted by the lasing mode alone (in this case: 
by an arbitrary but fixed mode) depicts a completely dif- 
ferent picture: The Fano factor diverges as the pumping 
is increased. This is easily understood by the qualita- 
tive description given above. For high pumping, every 
pump excitation eventually results in one photon being 
emitted from the cavity, but if there are several lasing 
modes the photon still has the freedom to chose one of 
those modes. These additional fluctuations can be that 
large that they eventually lead to a very large Fano fac- 
tor for large pumping. (It is obvious that the Langevin 
approach will break down eventually if the fluctuations 
become too large, as explained above.) 

The three test cases show that the model presented 
here is able to explain all relevant features of a laser. 



VII. RANDOM LASER 

A random laser is a laser where the feedback is not 
due to mirrors at the ends of the laser but due to chaotic 
scattering, either caused by scatterers placed at random 
positions or by a chaotic shape of the cavity jl], gj. If 
the mean outcoupling is weak a large number of modes 
in the cavity can be above threshold simultaneously [p5[ . 
As seen above, mode competition introduces additional 
noise into the modes. However, even if there are several 
modes above threshold, there only will be mode competi- 
tion if the modes are spatially overlapping and thus are 
"eating" from the same excitations. The main purpose 
of this paper is to answer the question whether in a ran- 
dom laser there is a relevant level of mode-competition 
noise or whether the radiation emitted in a laser line ap- 
proaches Poissonian statistics for strong pumping - both 
statements are mutually exclusively. 

We consider a chaotic cavity as depicted in Fig. |l| with 
a small opening to the outside. This problem becomes 
a stochastic problem by considering an ensemble of cavi- 
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ties with small variations in shape or scatterer positions. 
The coefficients appearing in Eq. (^) thus become ran- 
dom quantities. The statistics of these coefficients for a 
chaotic cavity with small opening is known |2^, ^] . The 
mean loss rate g of a cavity with volume V through a 
hole of diameter d at frequency u> is |3l| 



2-10"' 



2y2 



= 1-6 



(15) 



5 is the level spacing of the cavity. Its inverse 1/6 is 
the time needed to explore the entire phase space inside 
the cavity and can be identified with the round-trip time 
introduced for a "traditional laser" in Eq. (fi2|). 

In a chaotic cavity the modes ©i(r) can be modeled as 
random superpositions of plane waves |32| . This implies 
a Gaussian distribution for 0j(f) at any point r |3(|. The 
loss rate gi is proportional to the square of the gradient 
of 0(r) normal to the opening at the opening, hence its 
distribution is 



H9i) 



e -29i/g 



(16) 



and gi and gj are uncorrelated for i =/= j. 

It should be noted that the level spacing S is no random 
quantity, so that gi and U can be identified by choosing 
I /S as the unit of time. 

For simplicity we assume that the amplification pro- 
file w = 1 so that the distribution of the eigenfrequen- 
cies is not needed to compute Kij. (The distribution 
is known p9| so that an extension to non-constant w is 
straightforward.) 

Fig. |5| shows the computed Fano factor for a partic- 
ular sample from this ensemble (N p = 10, g = 0.5 but 
remember that the value of g of the lasing mode is much 
smaller than g Ell E3|). This kind of curve is typical 




FIG. 5: Fano factor of the radiation emitted from the (pri- 
mary) lasing mode (left axis, solid line) of some particular 
sample. The right axis (dashed line) depicts the number of 
modes above lasing threshold. Each additional mode crossing 
the threshold increases the Fano factor of the primary lasing 
mode. 
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FIG. 6: Fluctuations of the excitation density N(r) of the 
medium for the sample from Fig. ^. Depicted are the fluc- 
tuations (5 N 2 (r*i)) / N (f[) at the point f[ where the eigen- 
mode profile of the primary lasing mode has the largest 
magnitude (solid line, scaled by a factor 100), and the 
global quantities ([/ d 3 f 5N(r)] 2 }/ J d 3 rN(r) (long dashes) 
and (/ d 3 rSN 2 {r))/ J d 3 f~N{f) (short dashes). 



for all members of the ensemble, while the precise shape 
varies. When the first mode crosses the lasing thresh- 
old, the Fano factor goes through a maximum. While 
there is a global decrease with increasing pumping, addi- 
tional peaks are superimposed each time another mode 
crosses the lasing threshold. (In the following a mode is 
considered to be above lasing threshold if it contains at 
least 2 photons but the results are basically independent 
of whether one chooses 1, 2 or 10 photons.) The Fano 
factor approaches 1 plus some finite difference. Mode- 
competition noise thus gives a contribution to the noise 
but there still exists a lasing threshold that is well-defined 
by a peak of T . 

Similarly to computing the fluctuations of the Fano 
factor, it is possible to compute the fluctuations SN(r) 
of the density of excited atoms directly from Eq. <^). 
Fig. H depicts the computed fluctuations for the entire 
cavity (dashed lines) as well as for the point f\ where the 
eigenmode profile 0/(r) of the primary lasing mode has 
the largest magnitude. The former quantity peaks at a 
significantly larger pumping P which is immediately un- 
derstood by noticing that the primary lasing mode effects 
only part of the total cavity, and a significant part of the 
cavity is left "untouched" until more modes have crossed 
the lasing threshold. 

The two global quantities de- 

picted, ([/ d 3 fSN(f)} 2 )/ J d 3 rN(r) and 

(/ d 3 fSN 2 (f)) I J d 3 rN(f), differ by the inclusion 
of terms (N(ri)N(r2)), ?! ^ r-i- The different heights 
of the peaks (the first one is higher) demonstrate that 
(at least in the relevant interval of P, and on average) 
the density of excited atoms at different positions 
is positively correlated. This can be understood in 
the following simple picture: The photon densities 



7 




(^-l)/fl g/g Ni 

FIG. 7: The value of the Fano factor for the primary lasing mode depends on the number Ni of cavity modes above laser 
threshold, not on the other parameters. Unless otherwise noted, computed from ~ 9- 10° samples with g = 0.1. (a) Probability 
distribution of {JF — l)/g for g = 0.1, 0.2, . . . ,0.5. The five curves overlap almost perfectly, thereby demonstrating that the size 
of the opening does not influence the amount of mode-competition noise generated. (Computed from ~ 10 samples for each 
value of g with identical realizations for and gi/g for the five runs.) The inset shows the probability distribution from the 
large set with g — 0.1 plotted logarithmically, (b) Average of the Fano factor as a function of the outcoupling constant g of 
the lasing mode, (c) Average of the Fano factor as a function of N\. 



fii and the excitation densities N(r) are on average 
negatively correlated since each emission of an extra 
photon (Srii > 0) leads to the de-excitation of an atom 
[6N(r) < 0], and vice versa, hence (riiSN(r)) < 0. [This 
has also been confirmed by computing this correlator 
numerically from Eq. @]. Since the excited atoms at 
different positions communicate only via the radiation 
field, their density thus has to be positively correlated. 

It is difficult to relate the fluctuations of the exci- 
tation density of the medium to the properties of the 
emitted light. With increasing pumping, a peak of 
(5N 2 (ri))/N(n) starts to form (cf. Fig. §) at the same 
pumping that a peak starts to form for the Fano factor 
T (cf. Fig. ||) but the location of the maximum of the 
peak is significantly different for both curves. The com- 
plicated interplay between radiation modes and matter in 
a random laser does not allow for a simple understanding 
of the relation between these two quantities, and we will 
not discuss the fluctuations of the medium further in this 
paper since it focuses on the radiation properties. The 
complicated structure of the eigenmodes of a chaotic cav- 
ity is what makes a random laser fundamentally different 
from a "traditional" laser. 

In the following we will concentrate on the radiation 
and on the Fano factor far above threshold. P is chosen 
such that Pj ' g ss 10 7 (remember that the value of g of the 
lasing mode fluctuates). This is a compromise between 
a so large value as possible to ensure that the limiting 
value for P — > oo is approached as closely as possible, 
and a not too large value of P to avoid numerical prob- 
lems (remember that Fig. || already spans 11 orders of 
magnitude). 

The main results of a Monte-Carlo simulation with 
A^ p = 10 are depicted in Fig. 0. The scaled Fano fac- 
tor does not depend on the size of the opening (Fig. |7]a) , 
and only weakly on the outcoupling constant of the lasing 
mode (Fig. f7]b). As Fig. [?]c clearly shows, the true de- 
pendence is on the number N\ of modes above threshold. 
(The weak dependence of the Fano factor on the value of 



g of the lasing mode can be understood by noting that 
N\ is correlated with g of the lasing mode.) The finite 
value of T — 1 thus indeed is due to mode competition 
noise, as claimed above. 

For larger cavities, i. e. cavities with more modes in it, 
the distribution of ( T — l)/g changes from a peak near 
T = 1 to one that peaks at a finite value of (T — l)/g, 
as seen from Fig. As N p and N s increase, the effort 
to numerically compute the average solution from Eq. (j^) 
increases very fast, so that only a comparably small num- 
ber of realizations were computed (« 20000 for N p = 50 
and w 4000 for N p — 150), explaining the large sampling 
error in the histograms. [The speed could be increased 
significantly by developing an optimized algorithm for 
solving Eq. (^|).] For larger N p the average of (J 7 — l)/g 
becomes smaller as the large- T tail gradually disappears. 
(From N p — 10 to A^ p = 150 the average becomes smaller 
by about a factor 2; the average is difficult to compute 
since it sensitively depends on few samples with large J-.) 




0.2 0.4 0.6 0.8 1 
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FIG. 8: Distribution of the scaled Fano factor for cavities with 
different number N p — N s of modes. 
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FIG. 9: Small "virtual" cavities can be formed by scatterers 
in the random medium. Photons can be trapped very effi- 
ciently (i. e. small outcoupling) if the distances between the 
scatterers are compatible with the wavelength of the radia- 
tion) . 



VIII. INTERPRETATION OF EXPERIMENTS 

Experiments on random lasers are usually explained 
by the formation of small "virtual" cavities, which can 
"trap" laser light, so that it is scattered within a small 
volume many times before it can escape; see Fig. ^. (The 
linear dimension of such cavities was measured to be of 
the order of 100 wavelengths The chaotic cavity 

used as model in this paper should be understood as rep- 
resenting one of those virtual cavities. It is not obvious 
which values of the parameters (iV p , N s , gi, . . . ) are 
needed to explain the experiments. In the following we 
will argue that the important parameters are the average 
outcoupling g and, even more importantly, the probabil- 
ity distribution V(gi/g) as they together determine the 
number N\ of modes above lasing threshold. 

Above it was shown that N p and N s influence the Fano 
factor only weakly, i. e. only by a factor 2, and thus much 
smaller than the difference observed in the experiments. 
Even though it was not explicitly discussed in this pa- 
per, it is obvious that the choice of w(lo) and 0j(r) will 
not be important, either. This leaves g and V(gi/g) as 
parameters to explain the experiments. 

In this paper, a random laser is modeled by a chaotic 
cavity with a small opening. The size of the opening 
determines the average outcoupling <?, and all gi scale 
linear with g [see Eq. (|l6|)]. For a virtual cavity the 
average outcoupling cannot be computed in such a simple 
geometrical way. The outcoupling gi for the i-th mode in 
such a virtual cavity depends delicately on the positions 
of the scatterers and the wavelength of that mode. While 
no theory is available to compute gi or at least g for 
this case, it is likely that it will be relatively large as 
individual scatterers cannot be as effective as a massive 
wall with only one small opening. 

It was shown in Fig. 0a that T — 1 oc g. This is valid 
as long as the size of the opening is small compared to 
the square of the wavelength. If the opening becomes 
larger, the modes inside the cavity acquire a finite width 
(in frequency space) and start to overlap, severely com- 
plicating the theory p?| |, and it is not obvious how the 
behavior changes. Cao speculates that this overlap- 
ping prevents the formation of a fixed photon number in 



one particular mode as photons are constantly exchanged 
between modes with nearby frequencies. Furthermore, 
the Petermann factor of the lasing mode becomes signif- 
icantly large |22]| which might or might not increase the 
amount of fluctuations. While there is no proof that the 
amount of fluctuations is increased by these two effects, 
it seems to be obvious that the amount of fluctuations 
will not decrease due to them. Hence, T — 1 will at least 
increase proportional to the size of the opening — also 
for openings that are larger than the region of validity of 
the theory presented in this paper. 

The previous argument assumes that the number N\ of 
lasing modes inside a virtual cavity is the same as for a 
chaotic cavity with a small hole. Mode-overlap itself docs 
not change that number but for a larger opening the dis- 
tribution function V(gi/g) no longer has the form given 
by Eq. (fl6|). The form of V(gi/g) sensitively depends on 
the kind of outcoupling, and the number of lasing modes 
in turn sensitively depends on V(gi/g). For example, 
there already is a large difference between a cavity with 
one small hole and a cavity with two somewhat smaller 
holes (so that the total average loss rate is the same in 
both cases) p5| . It is very well possible that the form 
oiV{gi/g) may look significantly different from Eq. ( fl6| ) 
and could depend on many parameters of the sample. 

The differences in V{gijg) and thus in the number of 
lasing modes are thus the natural candidates to explain 
the differences observed in the two experiments. 

This prediction could in principle be checked exper- 
imentally by measuring the number of modes above 
threshold in one virtual cavity but to devise an experi- 
mental setup to do this seems very difficult, if at all pos- 
sible. The sample used by the group of Papazoglou || 
should have several spatially overlapping modes above 
lasing threshold (i. e. some modes above threshold are in 
the same virtual cavity), whereas in the sample by Cao j^] 
all modes above lasing threshold should be spatially sep- 
arated (i.e. be in different virtual cavities). One ex- 
planation could be that Cao's sample has more resonant 
feedback, so that the confinement of the lasing modes is 
stronger, compared to Papazoglou's sample. In the lat- 
ter, the modes would be extended over a much larger part 
of the sample (i. e. the virtual cavities are larger), giving 
them more possibility to overlap. 



IX. CONCLUSIONS 

In this paper we have developed a theory to compute 
the fluctuation properties of the radiation of a random 
laser. While for a standard single-mode laser the emit- 
ted radiation becomes coherent far above threshold, the 
radiation for a random laser is fluctuating more. It was 
shown that this extra noise is due to mode competition 
noise, i. e. due to the uncertainty of deciding into which 
mode to photon is emitted by induced emission. This 
noise is larger the higher the number of modes above 
lasing threshold is. 
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To be able to create mode competition noise, the com- 
peting modes have to be (at least partially) overlapping. 
On the other hand, if the profiles of the modes are over- 
lapping too much, usually only one of those modes will 
be above threshold. The amount of noise created thus is 
the result of a delicate interplay between these two com- 
peting effects. For a random laser modeled by a chaotic 
cavity filled with a laser dye, this leads to a finite in- 
crease of the Fano factor far above threshold, with the 
precise value depending on the number of modes within 
the cavity that are simultaneously above threshold for 
that particular realization of the disorder. In particular, 
the emitted radiation becomes coherent if only one mode 
is above threshold. 

Recent experiments on random lasers [)| gave con- 
flicting results on whether the noise is increased with 



respect to the Poissonian value. Even though it is not 
directly possible to model the differences in the two ex- 
periments, the theory presented in this paper suggests 
that this is due to the differences in the the number of 
modes above threshold. This number depends heavily on 
the specific system in question, so that the noise proper- 
ties of a random laser are not universal but depend on 
the (experimental) setup. 
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